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Effective stochastic equations for the continuous transitions of relativistic quantum fields in- 
evitably contain multiplicative noise. We examine the effect of such noise in a numerical simulation 
of a temperature quench in a 1+1 dimensional scalar theory. We look at out-of-equilibrium defect 
formation and compare our results with those of stochastic equations with purely additive noise. 
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I. INTRODUCTION 

Since phase transitions take place in a finite time, 
causality guarantees that correlation lengths remain fi- 
nite, even for continuous transitions. Because of the uni- 
versal presence of causality. Kibble and Zurek 0, Q 
suggested that it alone is sufficient to bound the size of 
correlated domains after the implementation of a contin- 
uous transition. There are several waysjsj of formulating 
causality bounds, but they all depend on the fact that, as 
the transition begins to be implemented, there is a max- 
imum speed at which the system can become ordered. 
For relativistic quantum field theory (QFT) this is the 
speed of light whereas, for superfluids, for example, it is 
the speed of second sound. 

The argument is very general. Consider a system with 
critical temperature Tc, cooled through that temperature 
so that, if T{t) is the temperature at time t, then r(0) — 
Tc- T'(O) = Tc/tq defines the quench time tq. Suppose 
that the adiabatic correlation length S,ad{t) = £.ad{T{t)) 
diverges near i = as 



t 



The fundamental length scale of the system is deter- 
mined from the microscopic dynamics. Although £,ad{i) 
diverges at i = this is not the case for the true 
non-equilibrium correlation length ^(t), which can only 
change so much in a finite time. Kibble and Zurek made 
the assumption that the correlation length ^ of the fields 
that characterizes the onset of order is the equilibrium 
correlation length f = ^ad(f) at some appropriate time t. 

For simple systems all estimates of t (the 'causal time') 
agree, up to numerical factors approximately unity 
Most simply, ^{t) cannot grow faster than c{t) — c{T{t)), 
where c{T) is the causal velocity at temperature T.[For 
relativistic theories c(T) — c, constant, whereas for con- 
densed matter systems we typically have critical slowing 
down, c(Tc) = 0.] This is true both before and after the 



transition. That is, i is defined by the condition that 
^ad{—t) ~ c{t) or ^ad{t) « " c(f). As a result, t is of the 
form 



t 



(1) 



where tq ^ tq is the cold relaxation time of the longest 
wavelength modes, and the critical exponent 7 depends 
upon the system. It follows that tq ^ t ^ tq. 

Domain formation, the frustration of the order pa- 
rameter fields, is often visible through topological de- 
fects, which mediate between different equivalent ground 
states. Since defects are, in principle, observable, they 
provide an excellent experimental tool for confirming this 
hypothesis when the possibility of producing them ex- 
ists. Kibble and Zurek made the further assumption 
that we can measure ^ experimentally by measuring the 
number of defects, assuming that the defect separation 
Cde/ = This identification of the initial domain size 

and defect separation then gives an estimate of the defect 
separation at formation of 



(2) 



where a = ^v. This is very large on the scale of cold 
defects which shrink to size £,ad{Tfin) = O(^o), where 
Tfin is the final temperature. We term a the Zurek- 
Kibble (ZK) characteristic index. 

In general, we find that the mean-field indices are 
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(3) 



for relativistic systems (e.g. weakly damped quantum 
fields), and 
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(4) 



for non-relativistic systems (e.g. superfiuids). There are 
exceptions to this rule, but we shall not consider them 
here. 
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II. STOCHASTIC EQUATIONS 

The question is whether these bounds, independent of 
the microscopic equations that govern the phase transi- 
tion, are remotely saturated in the physical world. In 
condensed matter physics several experiments have been 
performed to check ||2J). Although the results arc mixed 
the overall conclusion is positive ,4 i^j .fij i .Sj ■ 

Experiments cannot be performed for relativistic sys- 
tems and models of the early universe, when such tran- 
sitions were important, are too ambiguous to be helpful. 
In practice, an ab initio calculation from the microscopic 
field dynamics suggests the validity of the scaling laws, 
but does not permit an estimate of the efficiency of de- 
fect production outside the framework of mean-field (or 
large-N) approximations P.ITol|. 

In consequence, a more pragmatic check on the satu- 
ration of the Zurek-Kibblc bounds in relativistic systems 
has been numerical 11, 12], essentially from the empiri- 
cal damped relativistic Langcvin equation with additive 
noise, of the form 

^{x) + r]^{x) + {-W^ + m\t))(j>{x) + 2X(t)^{x)=^{x), (5) 

where, for simplicity, we consider the theory of a single 
real field. In Eq. the time-dependent mass ■m{t) — 
^~]^(t), that triggers the transition, is taken in the mean- 
field form 



m^{t) = -^,\{t), e{t) = il 



m 

Tr 



in the vicinity of t = 0, where is the cold mass which 
defines the symmetry-broken ground states, and e(t) ~ 
—t/TQ. For t > TQ it behaves as w?{t) 
equilibrium solution is 
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the symmetry breaking scale. Numerical simulations pTl 
IT^ show that the Kibble and Zurek scaling behaviour 
is recovered in the limits of small (underdamped) and 
large (overdampcd) rj respectively, although defects are 
produced with lower efficiency than anticipated in 

The overdamped case is the phenomenological non- 
relativistic time-dependent Ginzburg-Landau equation 



SF 

■q(j){x) = -— 



(6) 



for free energy _F, often appropriate for condensed matter 
systems. Equation Q arises from linear coupling to an 
environment. 

There is, however, a problem with justifying |SJ, with 
its assumption of linear dissipation, for underdamped 
relativistic QFT. Although linear couplings have played 
an important role in the history of decoherence and 
Langevin equations, additive noise alone cannot be jus- 
tified in QFT, where a pure linear coupling to the envi- 
ronment corresponds to an inappropriate diagonalisation 
of the fields. 



More generally, there are two mechanisms for inducing 
friction (dissipation) in a relativistic plasma: 

a) Changing the dispersion relation of existing parti- 
cles, as happens on varying rn?{t), leads to a change in 
scattering and decay rates, which leads to a change in the 
distribution of particle energies, which leads to friction. 

b) The creation of particles from the heat-bath, which 
leads to friction. 

In each case, we expect dissipative terms of the form 
^xj)^ to be important IT^ [l5l| . On the other side of 
the equation there is, equally, a problem with the noise 
in ||SJ) with regard to relativistic QFT. Noise is construed 
as a consequence of integrating out (tracing over) en- 
vironmental degrees of freedom, which are intrinsically 
non-linear. We assume Boltzmann statistics. In con- 
sequence, the noise, which guarantees that the system 
is ultimately driven to its ground states must, by the 
fiuctuation-dissipation theorem, contain a term of the 
multiplicative form 

This is confirmed in linear response theory which, for 
small deviations fro m equilibriu m, leads to Langevin 
equations of the form 



1 e quiii brmm, 
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9,,a^0(a;) - ^20-1- 2A03(x) + j d% Ki{x,y) c^{y) + 

+0(x) j d^y K2{x,y) 4>\y) = ^2 + <i>{x)S,i{x) + .(.Z) 

From the retarded nature of the Ks a more realistic, al- 
beit still idealised, equation than ((SJ for the real rela- 
tivistic scalar field (j) is 

0(a;) + (-V^ + m^{t))4){x) + 2A<?!)3(x) + 
+m<t>{x) + fii4>{x)4>'^{x) = ^i{x) + 0(x)e2(x), (8) 

where ^1,^2 are thermal noise. 

We note that, in ©, the multiplicative noise plays the 
role of a stochastic term in the temperature, on rewriting 
it as 

'(j){x) + (-V^ + ral{t))(j>{x) + 2\(j,^{x) + mH^) 
+T^2m<p\x)^Ux), (9) 



where 



in which 



T(t)-%(x). 



(10) 



For continuous transitions we are not aware of at- 
tempts to examine the effects of multiplicative noise for 
fields (as distinct from few degree-of-freedom quantum 
mechanics on scaling behaviour, 

even at the crudest phenomenological level. [This is not 
the case for discontinuous transitions, for which this pa- 
per is the counterpart to 26].] We see our work here as a 



3 



complementary study to the substantial analysis 0, ^| 
that has been undertaken for the simpler and often less 
believable equation to which it reduces for small 772 
and correspondingly weak noise. Finally, although we 
shall not pursue this here, we should not see ((SJ entirely 
in the context of QFT. There has been a considerable ef- 
fort in condensed matter to examine multiplicative noise 
in systems with simple ?70 dissipation, but with ^(f> multi- 
plicative noise (772 = C2 = 0) [23, llg ■ This leads to non- 
Boltzmann (Tsallis/Renyi) distributions [29l Isfiji. Mul- 
tiplicative noise can both induce phase transitions and 
restore symmetry . 



III. NUMERICAL SIMULATIONS 

It is sufficient for our purposes to restrict ourselves to 
a real field on the line (i.e. l-fl dimensions), for which 
= /dx^ . It was for such a system that the original 
equation (O was solved that case extension to 

higher dimensions gave few new insights as to the Kibble 
mechanism, and we expect the same here. 

Defects in this case are kinks. 



where m^(t) = -/z^e(t), /z^ = 1.0 and A = 1.0. and ^2 
are uncorrelated gaussian noise terms obeying 



fenkCa;) = iwtanh 



(11) 



of thickness ^0 — f^^^ and energy E = 0(/i^/A). Al- 
though, rigorously, there are no transitions for such short- 
range interactions in 1-1-1 dimensions, there is an effective 
transition for medium times. Typically, some time after 
the end of the quench, the field settles into a set of al- 
ternating positive/negative vacuum regions. These are 
separated by well defined kinks/ anti-kinks that evolve 
slowly, possibly annihilating each other for very long 
times. In this regime defects coincide with the zeroes 
of the field, making it straightforward to identify them 
in a numerical setting. Clearly this procedure is ambigu- 
ous for very early times, since zeroes occur at all scales, 
and only some of these will evolve into the cores of kinks. 
Here we will restrict ourselves to looking at kinks after 
the quench has terminated, thus avoiding any counting 
ambiguities. 



A. Numerical Setup 

To further clarify the possible effects of multiplicative 
noise terms in the mechanism of defect formation, we 
performed a numerical study of the model described in 
©. The approach followed is close to that in Ref.0 
for the case of a quenched one-dimensional system with 
purely additive white noise. We evolve the following 1-1-1 
Langevin equation: 



{Ux',t'Mx,t)) =2TSabSix' 



x)S{t' ~t), 



(13) 



aii^^i + 026 



+al](j)- m^{t)(j) - 2\f 



(12) 



The bath temperature T is set to a low value, typi- 
cally T = 0.01. The relative strength of the multiplica- 
tive and additive noise is measured by ai and a2, with 
corresponding dissipation terms obeying the fluctuation- 
dissipation relation. The values of ai and Q!2 vary be- 
tween different sets of runs, allowing us to compare their 
effects on the final defect density. 

When dealing with stochastic equations with multi- 
plicative noise one usually has to take into account that 
the continuous equation may not have a unique interpre- 
tation. This is the well known Ito-Stratanovich ambigu- 
ity |33| which is usually resolved by singling out a spe- 
cific discretisation of the equations of motion. It turns 
out that this problem has no relevance for the type of 
system we are considering. This could be seen by obtain- 
ing explicit Fokker-Plank equations for different time- 
discretisations of the model above. Since the multiplica- 
tive term depends only on the field (an not on its time 
derivative), the Fokker-Plank equation turns out to be 
the same for all alternative interpretations [2L ■23j . 

The equations of motion were discretised using a 
second-order leap-frog algorithm. We set 5x = 0.125 and 
(5t = 0.1 in a periodic simulation box with 8000 to 16000 
points. Note that the core of the defect (with finite tem- 
perature size l/rn) is resolved by 8 lattice points which 
should be enough for our purposes. 

Walls are identified by looking at zero crossings of the 
scalar field. As discussed above, this should be accurate 
for long times and low values of T - precisely the regime 
where we measure the final values for the defect density. 
For each individual quench, the final number of defects is 
determined by counting kinks at a final time, defined as 
a fixed multiple of the quench time-scale. There are more 
complex ways of defining a final defect density, namely 
by fitting the time-dependence of the kink number to 
exponential or power-law decay expressions. The sev- 
eral approaches were compared in The conclusion 
was that the straightforward kink counting performs less 
well in the very fast/slow quench-time limits, leading to 
slightly higher estimates of a. With this caveat in mind 
we will keep to the simpler approach since, as we will 
see, it is accurate enough to illustrate well the effects of 
multiplicative noise terms. 

For every fixed choice of ai and a2 we perform a series 
of quenches, with the quench times varying as tq =2", 
n = 1, 2, .., 9. Each quench is repeated several times (typ- 
ically 10) with different random number realizations, and 
the final defect number is averaged over this ensemble. 
The scaling exponent a can then be obtained by fitting 
the final defect number dependence on tq to a power-law 
of the form Atq"^ . 
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FIG. 1: Final defect density scaling power as a function of the 
additive noise strength. The multiplicative noise term is set 
to zero. Error bars represent the standard deviation of the 
result over 10 independent series of quench realizations. 



FIG. 2: Final defect density power versus additive noise 
strength. The multiplicative noise term is ai — 0, 0.5 and 
0.9 for the green (dahsdot), pink (dashed) and black (solid) 
curves respectively. Error bars as before. 



B. Simulation Results for Kink Densities 

We start with the simple purely additive noise case, 
corresponding to setting ai = in Eq. I|12|) . In FigQ] 
we can see the dependence of the scaling power cr, in 
terms of noise strength 012. Our results are very similar 
to those found in with a decreasing as the value of 
the dissipation, rj = of increases. This takes us from 
the relativistic regime where we expect cr ~ 1/3 to the 
overdamped case with a ~ 1/4. As observed by previ- 
ous authors, for very small values of the dissipation, the 
scaling fails to follow the power-law rule, a consequence 
of saturation. This explains the high, un-physical val- 
ues of a for a2 < 0.5 (corresponding to a dissipation of 
ry = 0.25). The quality of the power-law fit becomes poor 
in this parameter region, a further sign of deviations from 
the simple scaling behaviour. 

Next we look at how the introduction of multiplicative 
noise infiuences the results. In Fig|21we have the scaling 
power for a2 in the same region as FigQ] for three differ- 
ent values of multiplicative noise strength, ai = 0.0, 0.5 
and 0.9 respectively. Clearly the results change very lit- 
tle - within error bars, the three curves basically overlap 
each other. 

Though this result may seem disappointing at first, we 
should be aware that the significant range of variation 
of ai may differ considerably from that for a2- This is 
illustrated in Fig|21 where we have a for higher values 
of ai, with the contribution of the additive term set to 
zero. The pattern is the same as observed before for the 
case of additive noise. As ai increases, the scaling ex- 
ponent changes from relativistic values to those typical 
of an overdamped system. The transition between the 
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FIG. 3: Defect density power as a function of multiplicative 
noise strength ai. The additive noise is set to zero. Error 
bars as before. 



two types of behaviour takes place for values of ai con- 
siderably higher than a2- Defining the noise strength 
transition value as the one above which a < 0.3, we have 
a2 — 0.8 for purely additive noise and ai ~ 3.5 in the 
multiplicative case. 

These results can be understood if we make the simple 
assumption that the order of magnitude of the effective 
dissipation in the multiplicative noise case is given by 
PH - the mean of the term multiplying (j) in 
is the average of the square of the field dur- 



77 ~ a 
Eq.il 

ing the stage of the quench determining the scaling, that 
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haviour of the amplitude of the power-law. In Fig0] we 
show the value of log{A) where A is the amplitude of 
the fit for the final defect density Atq'^ . The results are 
shown as a function of the power a, corresponding to 
different ranges of ai and a2 in the purely multiplicative 
and additive noise cases respectively. As can be observed, 
A does not differ significantly between the two types of 
noise, for similar values of a. Physically this implies that 
not only the scaling power is similar in both cases, but 
also the overall amount of defects produced is not affected 
by the type of noise involved in the transition. Overall 
the conclusion seems to be that once we adjust for the 
value of the effective dissipation, the distribution of the 
defect number in the final configuration is independent 
of the properties of the noise terms driving the system. 
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IV. CONCLUSIONS 



FIG. 4: Amplitude of the power-law in terms of the power a 
for pure multiplicative (green dashed plot) and purely addi- 
tive (blue solid) noise. Error bars as before. 



is, slightly before the critical temperature is reached. If 
the typical value of {(jy^) is small, the effective dissipation 
for multiplicative noise should be reduced. As a conse- 
quence, relatively high values of ai should be required 
to take the system from the relativistic to the dissipative 
regime. 

This argument can be made more quantitative by not- 
ing that in the purely additive case the under-damped to 
over-damped transition takes place for rj ~ 0.8^ ~ 0.6. A 
similar effective dissipation with pure multiplicative noise 
would be reached for ai = 3.5 if {<j)^) ~ Q.e/S.S^ ~ 0.05. 
We measured the value of the mean squared field explic- 
itly in the simulations, and observed that at i = one 
has typically (0^) ~ 0.01 — 0.02, with the higher value 
corresponding to the slower quenches. This is indeed of 
the same order of magnitude as required by the above 
reasoning. We note that, if instead, we were to replace 
a\{(jP) by a^v"^ we would obtain an effective dissipation 
of = 6.1, one order of magnitude larger than the critical 
value. This suggests that even though kinks can only be 
identified with rigour once 



the period of the evolution responsible for setting the 
relevant defect separation scale takes place considerably 
earlier, in accord with (|T|l. 

The above result can be extended to systems with both 
kinds of noise, with the generic effective dissipation being 
given by 77 = 4- aKcjP). For ai < 1. the correction to 
the multiplicative component should be less than 0.05. 
This explains why the inclusion of multiplicative noise 
of this magnitude changes very little the additive noise 
result, as illustrated in Fig[21 

Finally we checked whether the inclusion of multiplica- 
tive noise leads to any appreciable change in the be- 



Since most effective equations of motion derived from 
field theory involve both additive and multiplicative 
noise, it is natural to wonder whether the presence of mul- 
tiplicative noise changes the dynamics of non-equilibrium 
phase transitions. Here we looked in detail at the forma- 
tion of defects in a 1-1-1 dimensional system undergoing a 
quench with both types of thermal noise present. Surpris- 
ingly, we found that the properties of the defect popula- 
tion after the quench can be well described in terms of the 
Kibble-Zurek scenario, if we take into account the effects 
of multiplicative noise in the dissipation. In particular, 
multiplicative noise terms increase the dissipation by an 
amount of the order of the noise amplitude, times the 
value of the mean field square at the transition time. As 
in the purely additive case, we observe that for low values 
of the effective dissipation the system behaves in a rela- 
tivistic fashion, with the final number of defects scaling 
as a power of 1/3 of the quench time. For higher values 
the system enters an over-damped regime characterized 
by a lower scaling power, nearer to 1/4. Although our 
multiplicative noise was the most simple (tracing over 
short wavelength modes will also induce (lP£,3 noise) this 
suggests that the Kibble-Zurek scaling laws are robust. 

It is interesting to note that these results still leave 
open the question of whether the defect density (or equiv- 
alently, the freeze-out correlation length) is set before or 
after the transition takes place, i.e. for t < or i > 0. 
Strictly, the value for the average field square ((/)^) used in 
defining the effective dissipation, should be evaluated at 
t. Unfortunately the value of this quantity varies slowly 
in the vicinity of the critical point and as a consequence 
it is not possible to determine whether the correct value 
is fixed above or below Tc. We can nevertheless, use our 
model as a setting for answering this question, by per- 
forming quenches where the value of the effective dissipa- 
tion is forced to change att = 0. By looking at the defect 
scaling in cases where the shift in the dissipation takes 
the system from an under-damped to an over-damped 
regime at t = 0, the significant stage of the evolution 
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should become apparent. We will clarify these points in 
detail in a future publication. 
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